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For the first time, the cumulant 2-body reduced density matrix (= 2-matrix) of the 
spin-unpolarized homogeneous electron gas (HEG) is considered. This 7c proves 
to be the common source for both the momentum distribution n{k) and the static 
structure factor S{q). Within many-body perturbation theory, this 7c is given by 
only linked diagrams (with 2 open particle-hole lines as well as with closed loops 
and interaction lines). Here it is worked out in detail, how the 1-body quantity n{k) 
follows from the 2-body quantity 7c - through a certain contraction procedure, cf 



Eqs.( 2.25 )-( 2.28). In particular, this 7c is developed for the high-density HEG. Its 
correctness is checked by deriving from it n{k) and S{q), known from the random- 
phase approximation (RPA) . This study opens the way to a more sophisticated HEG 
description in terms of cumulant geminals or/and variational methods. Besides, the 
cumulant structure factor (CSF) of the exchange in lowest order is explicitly given 
and sum rules for the CSFs and their small- and large-g behavior (beyond RPA) are 
systematically summarized within the plasmon sum rule, coalescing theorems, and 
the inflexion-point trajectory. 

List of Symbols/Shorthands/Abbreviations 

RDM reduced density matrix, SR sum rule, p homogeneous electron density 

MD momentum distribution n(k), f{r) 1-body RDM 

PD pair density g{r), CPD cumulant PD h{r): g{r) = 1 — \f'^{r) — h{r) 

SF structure factor S{q), CSF cumulant SF C{q): S{q) = 1 - lF{q) - C{q) 

"a" antiparallel spin: ga,{r) o S'a(g), Kir) O Ca(g) 

"p" parallel spin: gp{r) o Sp{q), hp{r) ^ Cp{q) 

FT Fourier transform, F{q) = FT of f'^{r), FS Fermi surface |fc| = 1 
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I. INTRODUCTION 

The lowest-level quantum-kinematics of an extended many-electron system (within the 
Born-Oppenheimer approximation) is contained in 1- and 2-body densities and reduced 
density matrices (RDMs): pi(l)= electron density, 7i(l|l')= 1-body RDM (1-matrix), 
P2(l, 2)= pair density (PD), 72(1 2|2')= 2-matrix. The well-tried and widely used density 
functional theory (DFT) yields pi (but not 71 and P2). If a density-matrix functional theory 
would exist, it would yield 71 (but not P2); if a pair-density functional theory would exist, 
it would yield p2 (but not 71). If an effective 2-body scheme would exist, which yields 
approximately the 2-matrix 72, then from it would follow both the PD p2 (by taking the 
diagonal elements) and the 1-matrix 71 (by means of the contraction procedure 2' = 2 and 
J d2 72(11 1', 2|2)). For extended systems this requires the cumulant decomposition of 72. In 
this paper (it continues from [5lll2l|53] and is related to [361138] ). the cumulant 2-matrix 
7c (dimensionless: x) is developed as a decisive key quantity and its properties are studied 
for a "simple" model system, namely the high-density spin-unpolarized homogeneous 
electron gas (HEG). 

Although not present in the Periodic Table, this HEG is still an important and so far 
unsolved model system for the electronic structure theory, cf e.g. [Ij. Its advantage: pure 
"correlation" and no "multiple-scattering". In its spin-unpolarized version, the HEG ground 
state (GS) is characterized by only one parameter r^, such that a sphere with this Wigner- 
Seitz radius contains on average one electron [2]. It determines the Fermi wave number 
as kp = l/{ars) in atomic units (a.u. ) with a = [4/(97r)]i/3 ^ 0.521062 and it measures 
simultaneously both the interaction strength and the homogeneous electron density p, such 
that high density corresponds to weak interaction and hence weak correlation [3]. For recent 
papers on this limit cf |1HT2]. = corresponds to the ideal Fermi gas. Its GS energy per 
particle is Cq = 3/10, measured in units of kp. The Coulomb repulsion v{q) = q^^jq^ with 
q1 = Aars/n causes deviations. This is in lowest order the exchange energy Cx = — (3/47r)ars. 
But in next order the long range of the Coulomb repulsion makes the correlation energy Ccorr, 
defined by e = eo + Cx + Ccorr + 0{rs), to behave non-analytically at the high-density limit: 
Ccorr ~ r^[lnrs + const + 0{rs)] p^HTSl [38 j . This non-analytical behavior of e carries over 
to the kinetic and potential components, t respectively v, through the virial theorem after 
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March (1958) [IS]: 



f = r^— e , t = -r- e. (1.1) 

dr ., aVs Tg 



The energy components t and v follow from the simplest quantum-kinematical quantities, 
namely n{k) and S{q), the momentum distribution (MD) and the (static) structure factor 
(SF), respectively: 



oo oo oo 



t = I d{k') n{k)— , I d{k') n{k) = I , I d{k') n{k)[l - n{k)] = c , (1.2) 





oo 

V = — 





oo oo 

y"f^[l-5(g)]^, 5(0) = 0, ld{q')[l-S{q)] = l + c,. (1.3) 



k and q are measured in units of kp and r in units of k^^. The quantity c is hereforth 
referred to as Lowdin parameter, because Lowdin was the first one who queried the meaning 
of the trace of the squared 1-matrix. Because of < n{k) < 1, c measures the r^-dependent 
non-idempotency of n{k), being zero for = (ideal Fermi gas) and increasing with 
Tg. Thus c measures simultaneously the strength of correlation. Note its particle-hole 
symmetry. Ci is another (short-range) correlation parameter, also vanishing for — )■ 0. 

Describing electron pairs, one has to distinguish SFs S'a,p(q') for pairs with antiparallel spin 
("a") and with parallel spins ("p") and corresponding pair densities (PDs) 5'a,p(^); where 
ga{r) and gp{r) describe the Coulomb hole and the Fermi hole, respectively. The 5'a,p('") and 
Sa.,p{q) are mutually related through Fourier transform (FT): 



^a(r) - 1 



oo oo 

[diq'f^SM ^ = a' [ d{r'f^ lig^ir) - 1] , (1.4) 

J qr J qr 2 



oo oo 

,p(r) - 1 = / d{q'f^ [Sp{q) - 1] O 5p(g) - 1 = a' [ d{r' f^ \[gp{r) - 1] . 
J qr J qr z 



From this follow the "total" PD g{r) = [g^{r) + ^p(r)]/2 and the "total" SF S{q) = S^{q) + 
Sp{q): 

oo oo 

g{r)-l= [d{q'f-^l[Siq)-l] ^ S{q) - 1 = a' [ d{r'f-^ [ g{r) - 1] . (1.5) 
J qr z J qf 



5'a,p(0) = fix the normalizations of 1 — 5'a,p('"), known as the perfect screening SR, whereas 
the coalescing (or on-top) values ga.{Q) and 5'p(0) = fix the normalizations of S'a(g) and 
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Sp{q) — 1, respectively. ga{0) = 1 — ci or g{0) = (1 — ci)/2 shows the physical meaning 
of Ci- It determines the short-range correlation in terms of 5'a(0), the on-top value of the 



Coulomb hole ga{r) with < ga{0) < 1. For further details (normalizations) cf (B.l), (B.2). 



In view of Eqs.(1.2) and (1.3), one may ask which pecuharities of n{k) and Sa,p{q) cause 
the non-analyticities of t and v, respectively. As shown in 0, the above-mentioned 
drastic changes, when switching on the long-range Coulomb interaction, show up in the 
redistribution of the non-interacting MD nQ{k) = G(l — k) within thin layers (thickness Qc) 
inside and outside the Fermi surface |fc| = 1 and with a remaining finite jump discontinuity 
at |fc| = 1: zp = n{l~) — n{l^), < Zp < 1. They show up also in the behavior of 
S{q) within a small spherical region (radius Qc) around the origin of the reciprocal space 
with the plasmon SR S{q ^ qc) = g^/(2a;pi) [T71 [18], which causes an inflexion 

point ginfl, •S'infl ~ Wpi, where Upi = qc/\/3 is the plasma frequency, measured in units of 
kp. The mentioned non-analyticities and redistributions let ordinary perturbation theory 
fail, e.g. with S{q — )■ 0) ~ 1/q and n{k) = noik) + An{k) with nQ{k) = 0(1 — k) and 
An{k — )■ 1^) ~ ±l/(^ — 1)^- That the GS energy e{rs) diverges in 2nd order, has been 
shown by Heisenberg (1947) [13]. Macke (1950) [14j has repaired this failure by means of 
an appropriate, physically plausible partial summation of higher-order perturbation 
terms for the GS energy per particle e{rs), a way which has been developed further by 
Gell-Mann and Brueckner (1957) [T3]. This so-called ring-diagram summation, cf Fig.l, is 
also known as the random phase approximation (RPA) . This summation has been developed 
for n{k) by Daniel/Vosko (1960) [I9l and Kulik (1961) [20], an analytical extrapolation 
for n{k) is given in [211 122] , for the spin-polarized case see [23], recent quantum Monte 
Carlo calculations of n{k) for = 1, ■ ■ ■ , 10 are in [23], zp for high-densities is in [25] . 
zp calculated for < 55 is in [26j. The ring-diagram summation for S{q) has been done 
by Glick/Ferrell (1960) [27], Geldart (1967) [2B], and Kimball (1976) |29j. For attempts 
to go beyond RPA cf eg [231 ED]- The on-top behavior of g^^p (describing the short-range 
correlation) influence the large-g behavior of Sa,p{q) and n{k). This comes from theo- 
rems referred here as the coalescing cusp and curvature theorems, cf [31] and [32] and App.B. 



With n{k) also t is available (but not v) and with S{q) it is available v (but not t). Does 
a quantity exist which contains both n{k) and ^(g)? "Old" and recent attempts in that 
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direction (with variational or direct calculations of the 2-matrix) are in [59l [65| 166] . In [12] 
it has been shown that the self-energy S(/c, u) as a. functional of t{k) = and v{q) = q^/q^ 
is such a quantity [69j. Here it is shown - using the concept of reduced density matrices 
(RDMs) - that n{k) and S{q) have their common origin in a quantity called 'cumulant 2- 
matrix', x(l|l',2|2') with the symbolic variable 1 = (ri,o"i). This matrix appears if one 
tries to represent the 2-matrix of an interacting system in terms of the 1-matrix (according 
to an independent-particle or Hartree-Fock model). There remains an irreducible part, 
which can not be reduced due to correlation. For the PD and for the SF, these cumulant 
decompositions are 

g.{r) = 1 - K{r) , = -CM , (1-6) 

g^{r) = l-f{r)-h^{r), S,{q) = 1 - ^F{q) - C,{q) . (1.7) 

Note the asymmetry with respect to "a" (Coulomb hole) and "p" (Fermi hole). With the 
"total" cumulants h{r) = [h^ + ^p(^)]/2 and C{q) = Ca(g) + Cp{q) it follows 

g{r) = 1 - ^f{r) - h{r) , S{q) = 1 - ^-F{q) - C{q) . (1.8) 

/(r) is the dimensionless 1-matrix [following from n{k), cf (2.4)] and F{q) - referred to as 



HF-function - is the FT of P{r), cf (B.15). hs_^p{r) are the cumulant PDs (CPDs) and their 



FTs are the cumulant SFs (CSFs) Ca,p(g): 

oo oo 

K{r) = ! d{q'f^ CM ^ CM = [ d{r'f^ IhM , (1-9) 
J qr J qr 2 



oo oo 

hM = I d{q'f^ CM ^ CM = 0? j rf(r^)^ \K{r) ■ (I-IO) 



Within RDM theory one can show that the cumulant 2-matrix x(l|l') 2|2') is the source for: 
(i) the CSF C{q) [or equivalently the CPD /i(r)] and (ii) also for the MD n{k). Whereas 

(i) results simply from x(l|l,2|2), i.e. taking the diagonal elements V = 1, 2' = 2, the case 

(ii) follows from the slightly more complicated RDM contraction of x(l|l') 2|2') with 2' = 2, 



/ d2 x(l|l',2|2) and FT. The Eqs.(1.9) and (1.10) contain the CSF properties (B.l) and 



(B.2) needed below. 
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General remarks on RDMs and cumulants: Whereas finite many-body systems can 
be described by many-body wave functions $(1,2,---) as solutions of a Schrodinger 
equation, extended systems have to be described by a hierarchy of A^-representable RDMs 
as solutions of the BBGKY-like hierarchy of contracted Schrodinger equations. And: 
whereas for finite systems the concept of cumulant (2-body, 3-body, ...) matrices can be 
applied, for extended systems it is a must, because - unlike the RDMs - all these cumulant 
matrices are size- extensive entities, what allows the thermodynamic limit with — > oo, 
Q oo, p = N/Q = const. Closely related with this property is that they are given within 
perturbation theory by non-vacuum linked Feynman diagrams. In App.A, the systematic 
definition of cumulant matrices in terms of generating functionals is summarized. Recent 
(quantum chemical) papers on RDMs, cumulants, contracted Schrodinger equations, 
correlation strength, correlation entropy, quantum entanglement. Berry phases etc. are 
[55]- [6l] and refs. therein. [65] and [66] deal (for finite systems) with variational calculations 
of the 2-matrix. Recent HEG papers are [HZ] and [HE] and refs. within. 



Here is roughly sketched, what will be presented in Sees. II, III, IV in more detail. In 
this paper it is aimed to present x the RPA with all those terms of x contributing 
to the correlation energy Ccorr terms up to Tglnr^ and r^. The correctness of this x is 
controlled/checked by deriving n{k) and S{q). A comparison is performed between these 
'new' RPA results with the 'old' ones of Daniel/Vosko [IH], Kulik [20], and Kimball [22], 
respectively. For this purpose the diagrams of Figs. 1-5 are needed. Fig.l shows the 
Yukawa-like screening of the bare (long-range) Coulomb repulsion v{q) = ql/q^., replacing 
it by the frequency-dependent interaction v{q,rj) = v{q)/[l + v{q)Q{q,ri)], where Q{q,ri) is 



the particle-hole propagator (C.3). The bare Coulomb repulsion shows up in divergencies: 
for the correlation energy Ccorr see [l3], for n{k) and S{q) see [6]. These divergencies are 
eliminated through Macke's partial summation of diagrams. Fig.2a shows the interaction 
of 2 particle-hole lines running from 1' to 1 and from 2' to 2. This is called "d" = direct 
diagram to distinguish it from "x" = exchange diagram of Fig. 3a with one line running 
from 1' to 2 and another one from 2' to 1. These lowest-order RPA terms of Figs. 2a 
and 3a are called Xdr and Xxr, respectively, "r" = ring-diagram summation of Fig. 1. The 
cumulant matrices Xdr,xr yield the CPDs /^dr.xr and the CSFs Cdr.xr of Figs. 2b and 3b and 
the interaction energies of Figs. 2c and 3c. But how to obtain An{k)7 Using the self-energy 
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oj), the lowest-order S - diagrams are in Figs. 4c and 5c: An{k) = fij-lk) + n^{k) + ■ • • . 
Of course, RDM-theory must end up with the same results. Indeed, nj.{k) of Fig.4c comes 
from Xxr of Fig.Sa by the contraction 2' = 2 and J d2, indicated in Fig.4b by a small 
circle with the final result rij-lk) of Fig. 4c. Similarly, nx{k) of Fig.5c comes from the ladder 
diagram (in its exchange version) of Fig.Sa, contraction makes Fig. 5b with the final result 
nx{k) of Fig. 5c. Contraction in terms of diagrams mean to transform a 2-line diagram into 
a 1-line-diagram. This study opens the door for variational 2-matrix calculations and for 
discussing HEG in terms of cumulant geminals and their weights. 



The outline of the paper is as follows. In Sec. II the basic definitions of RDMs, their 
cumulant decompositions and the contraction SR are presented. It is derived in detail how 
the 2-body quantities S{q) and g{r) with its cumulant pedants C{q) and h{r), together 
with the 1-body quantity n{k), follow from the cumulant 2-body matrix x- The decisive 



contraction procedure is in Eqs.(2.25)-(2.28). In Sees. Ill and IV the lowest-order cumulant 
2-matrices and what follows from them are presented. Sec.V gives a summary and an out- 
look. App.A explains the exponential-linked-diagram theorem as the base for a systematic 
cumulant decompositions. App.B deals with SF properties in terms of the coalescing cusp 
and curvature theorems, of the plasmon SR, and of the inflexion-point trajectory. Table 1 
and 2 contain the small- and large g-behavior of the CSFs and SFs (similar as in [36] )• In 
App.C the particle-hole propagator, the Macke function I{q) and contour integrations are 
described. App.D lists characteristic correlation parameters, vanishing for — )■ 0. 



II. BASIC DEFINITIONS AND RELATIONS 

2.1 The 1-matrix and the 2-matrix: definitions and properties 

The starting point is the definition of the 1-matrix 71 and of the 2-matrix 72 for the HEG-GS 
in terms of creation and annihilation operators, ipi and ipiy respectively: 

7i(l|l') = , 72(1|1',2|2') = (^1,4^2^1) . (2.1) 

The shorthands 1 = (ri,(Ti) and j = ^ J d^^i are used. The hermiticity is obvious and 
Fermi symmetry means 72 — ?■ —72 if 1 and 2 or 1' and 2' are interchanged. N = J dl iplipi 
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is the total particle-number operator. The bra vector " (" as well as the ket vector ")" are 
A^-body states, therefore ipi) and {ipl are (A^ — l)-body states, thus 

dl = N , j dld2 {^li;l^2^i) = N{N - 1) (2.2) 

holds for the normalization. The contraction SR 

d2{44^P2i^r)=liil\l'){N-l) (2.3) 

describes how the 1-matrix 71 for finite N results from the 2-matrix 72. 

It follows the spin structure of 71, the dimensionless 1-matrix /(r) and its FT: 

7i(l|l') = ^6^„a'Mri\r[) , 7i(^ik'i) = pf{rir) , ni' = \ri - r[\ , (2.4) 

f{r) = J d{k') ^ n{k) , /(O) = I d{k^) n{k) = 1 , 


/'(O) = 0, /"(O) = -2t/3 (with to = 3/10 and /^'(O) = -1/5 for r, = 0). Its natu- 
ral orbitals are plane waves, because the system is extended and homogeneous. n{k) is 
the momentum distribution (MD), k and r are measured in units of kp and /cp^, respectively. 

The symmetric and antisymmetric 2-body spin functions 



S±{ai\(T[,a2\cr'2) = ^((5^j,^^(5^2,^^ ±5^j,^^5^2_^/) with ^ ^ 5±(cri|(Ti, cralcxa) = | 



enter the spin structure of 72(1|1', 2|2') as 

72(1|1', 2|2') = 5_(o-i|ai,a2|cr2) ltiri\r[,r2\r'2) + 5+{ai\a[,a2\(T2) I2 ifi\r[,r2\r'2) (2.6) 



[231 IM] . It defines - together with the basic definition (2.1) - the singlet 2-matrix 7^ and 
the triplet 2-matrix 7^, being (in space) symmetric and antisymmetric for the interchanges 
Vi ^ 7^2 or r[ ^ r'21 respectively. The normalizations 

f^'rirf\2 72^(^11^1,^-2^2) = %i%±\\ (2.7) 



2 V 2 

follow from the contractions 
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which are in their turn consequences of the total normahzation (2.2) and the singlet-triplet 



representation (2.6). Eqs.(2.7) and (2.8) show the necessity to find an equivalent writing 



for 72, which allows the thermodynamic limit with — t- oo, Q — t- oo, p = N/Q = const. 
This is just done by means of an additive decomposition of 72 in a Hartee-Fock like term 
7hf, which is defined in terms of the (full, exact) 1-matrix 71 only, and a non-reducible 
remainder - the cumulant 2-matrix 7c. Thus 72 = 7hf — 7c as the first step of a more 
general expansion, see App.A. 

2.2 The HF approximation 

The 2-matrix in a Hartree-Fock like approximation (in terms of the exact 1-matrix only) is 



7hf(1|1',2|2') =7i(l|l')7i(2|2')-7i(l|2')7i(2|l') 



(2.9) 



This may by also addressed "exact exchange". Using Eq.(2.4), the same structure (2.6) as 
for the total 2-matrix 72(1|1', 2|2') results: 



7hf(1|1',2|2') = 5^{ai\a[,a2\a2h^piri\r[,r2\r2) + 5+(o-i|o-i, o-2k2)7HF(^il^i' ^2^2) 



(2.10) 



with 



7HF(^lkl>^2|T'2) = 71 (^lk'l)7l(^2|T'2) ±71 (^1^2)71 (^2|t'i) 

= ^ [/(rii')/(r22') ± /(ri2')/(r2i')] ■ 



Their normalizations 
contain the above introduced Lowdin parameter c, which is originally defined by 



2 V 2 2 



(2.11) 



(2.12) 



l-c=| j dh,dh2\\li{ri\r2)? = a' j d{rl^)\f{r,2) = j di^) n\k) . (2.13) 



Note p J d^Ti = pQ = N and note that the homogeneous density p, measured in units of 
kl, equals 1/{37t^) and that d^r^/iSTf^) = a^dirf^), a = [4/(97r)]i/l 



The contractions of 7^p are 



:7il^i|^i. 



\ 1 



(2,14) 
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where 



c(ri r 



P 



4:71/3 



c{k) e 



, c{k) = n{k)[l - n{k)] , / rfV c{r\r) = Nc . (2.15) 



c{k) and c{r,r') are referred to as "Lowdin function" and "Lowdin matrix", respectively. 
[c{k) appears in [62J as — 2efc]. 



2.3 The cumulant decomposition and cumulant geminals 

With this prelude, the above mentioned decomposition 

72(1|1',2|2') =7hf(1|1',2|2')-7c(1|1',2|2') 



(2.16) 



defines the cumulant 2-matrix 7c. Its spin-structure is similar as in (2.6) and (2.10) 



7c(l|l',2|2') = 5-(^ik;,^2k^) iUr,\r[,r2\r',) + 6+{a^\a[,a2W',) % {n\r[,r,\r',) (2.17) 
and it defines (through 7^ = 7^p — 7^ ) the cumulant 2-matrices 7^ with the normalization 

d?rid^r2 7^ (^i^i, ^2^2) = =F^c , (2.18) 



following from the normalizations (2.7) and (2.12) . This ensures the normalization of 7, 



to be / dld2 7c(l|l,2|2) = Nc, as it should. With the help of Eqs.(Q and \2.1A\ the 
contractions of 7^ are 



j d^r2 7^(^ik'i,^2k2) = =F^ c{ri\r[) . 



(2.19) 



Eqs.(2.18) and (2.19) show that the Xc size-extensively normalized and contracted. With 



7^ = -^X±: the dimensionless cumulant 2-matrices x± are introduced. Their normalization 



follows from (2.18) as 



00 

J dirl^) X±{ri\ri, rsjra) = t2c . (2.20) 



The introduction of the singlet /triplet matrices 7^ and x±i is only a necessary intermediate 
step to present its spin-parallel/antiparallel components Xp,a in terms of linked diagrams. 
For diagonal spins ci^ 2 = '^1,2 it is 5±{ai\ai, cr2|(T2) = (1 ± 5o-i,<t2)/2- So for antiparallel spins 
(cTi = —0-2, hence 5o-i,o-2 = 0), respectively parallel spins (ai = (72, hence 5o-i,cr2 = 1) the spin 



structure (2.17) leads to 



(2.21) 
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The arguments of these matrices are {ri\r[,r2\r2). Their diagonals {ri\ri,r2\r2) determine 
the normahzations 

oo oo oo 

j dirl^) ^Xa = 0, J dirf^) ^Xp = c ^ j d^rl^) ^x± = Tc . (2.22) 



The singlet /triplet matrices x± define symmetric and antisymmetric cumulant geminals 
V'J(ri,r2) with corresponding weights z/^ and with T^ri^k — "Fc as total weights. Which 
properties do these geminals have? Are they discrete (bound) or/and continuous (scatter- 
ing) states, the latter with phase shifts? Is there a kind of an aufbau principle for the z/^? 



2.4 Linked Diagrams 

As mentioned above, the cumulant 2-matrices x± cire given by linked diagrams. To each 
symmetrized direct diagram di=Xdi (^i|^i, '^21^2), d2, ■ ■ ■ with 2 open particle-hole lines 
(one running from r[ to Vi and another one from r'g to r2) belongs an exchange diagram 
Xi=Xxi(^i|^'i7 '^21^2) — Xdi (^i|^27 '^21^1), X2, ■ ■ ■ , such that for r'^ = r2 and / d^r2 there is 
only one line running from r[ to Vi (like the diagrams of the self-energy S, cf [12j). With 
these building elements Xd = Xdi + Xd2 + ■ ■ ■ and Xx = X^i + Xx2 + ■ ■ ■ the singlet / triplet 



components are X± = Xd ± Xx and from (2.21) follow the decisive relations 



Xa 



Xd 



Xp = Xd - Xx r\ X± = Xd ± Xx 



(2.23) 



again with the arguments {ri\r[,r2\r2)- The spin-averaged sum x can be written as 



X=;^[Xa + Xp] or x 



The a- and p-components are equally weighted, whereas the singlet- and triplet-components 
have the weights 1/4 and 3/4, and the d-and x-components have the "weights" 1 and -1/2, 



Xd 



:Xx 



(2.24) 



respectively. From Eqs.(2.23) and (2.24) follow corresponding relations for the PDs g.. 



a,p 



and their cumulant correspondings /ia,p- Note that for their FTs S'a,p and Ca,p it holds 
5 = 5a + 5p and C = Ca + Cp. 



Now the question is: how real physical properties like the PD's 5'a,p(^) (by taking the 
diagonal elements) as well as the MD n{k) (by Fourier transforming the contracted 
off-diagonal elements) follow from the (more formal) key quantities Xa,p, respectively Xd,x- 
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2.5 Pair Densities and Structure Factors 

With the cumulant PDs /ia('"i2) = Xa(^'i|^'i5 ''"21^2) and hp{ri2) 



Xp[ri\rur2\r2} 



and 



their Fourier transforms Ca^p{q) following from (1.9) and (1.10), the PDs ga,p{f) and 



SFs 5'a,p(g) are given by (1.6)-(1.8). This shows that for the spin-parallel quantities, the 
functions /^(r) and its FT F{q) are needed, cf App.B. So from the cumulant 2-matrices 
Xa,p(^i|^'i, ''"2|^"2) follow not only the spin-antiparallel PD ga{r), but also the spin-parallel 

PD gp{r) (supposed the MD n{k) has been calculated first, cf next paragraph). Rigorous 

00 

theorems for Ca,p(g — )■ 0), Ca,p(g — 00), and / d{(l^)Ci,^p{q) are summarized in App.B. 



2.6 Contraction Sum Rule and Momentum Distribution 

In terms of the Lowdin function c{k) = n{k)[l — n{k)] and of the contracted cumulant matrix 
(the notation means r2 = r2 and J d^r2 ■ ■ ■ ) 



d?r2 
2-37r2 



X±{.ri\r\,r2\r2) r\ x±ik) 



2-37r2 



X±(^i|^'i) 

(2.25) 



the contraction SR (2.19) can be written as x±{k) = ^c{k). If the expressions (2.21) and 



(2.23) are contracted and Fourier transformed, then 
1 



Xdik) 



-c{k) + c{k)] = , Uk) 



1 



This is the contraction SR, which allows to calculate n{k) from 

d^rw -ikr,,, f d^r2 



c{k) - c{k)] = -c{k) . 



(2.26) 



2-371^ 



XA'ri\r„r2\r2) 



(2.27) 



supposed Xx is approximately known. The 'direct' diagrams Xd(^) with 2 lines do not 
contribute to n{k), which alone derives from the 1-line diagram x^{k). If no{k) = 0(1 — k) 
means the MD of the ideal Fermi gas, then n{k) = nQ{k) + An{k) defines its interaction 
induced correction An{k) ^ for k ^ 1 being normalized as J d^k An{k) = 0. Thus the 



contraction SR (2.26) reads finally 



T An{k) + [An{k)f = x^{k) for k^l 



(2.28) 



what simplifies as An{k ^ 1) ~ TXx(^) small correlation tails An{k), as this is the 
case for the RPA with Xx{k) = Xxr{k) + ■ ■ ■ and An{k) = Arij-lk) + ■ ■ ■ . The ellipsis mean 
beyond- RPA terms. From a diagrammatical point of view one may expect, that already 
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is the MD, but this is not the case: as derived above, a factor sign(l — k) has to be 
added, to transform Xx(^) to An{k). So contraction means r'g = r2, J (Pr2 and adding a 
factor sign(l — k). - From An{k) follows the Lowdin parameter c, which fixes some CSFs 
at long wavelengths: Cx(0) = — c, Cp(0) = +c, and C(0) = c. 



2.7 Energies 

oo 

From An{k) follows the kinetic correlation energy At = J d{k^)An{k)k'^/2. The same An{k) 



makes AF = F—Fq with (B.18). This allows to calculate the interaction energy v = vuf+Vc, 
which is built up from the HF part whf = + Avuf and the cumulant part Vc- 



3 ■ 4 2 



Fo{q)v{q) , AvuF = - 



d{q^) 1 
3-42 



AF{q)v{q) , Vc 



d{q') 
3-4 



C{q)v{q) . 
(2.29) 

— ars/(47r/3) is the HF energy in lowest order. C = Ca+Cp is the FT of h = [hs,+hp]/2, 



cf (1.9) and (1.10). The cumulant interaction energy Vc = Vs_ + Vp is part of the interaction 
correlation energy fcorr = ^?^hf + ^^c- Thus, the energy-components t and v are available, 
defining the total energy e = t + f as a function of with v = rgde/drg and t = e — v. The 
correlation energy is Ccorr = At + Avuf + Vc- 



2.8 Summary of Sec. II 

If X±(^i|^i) '^21^2) is given in terms of diagrams di,d2,--- with 2 not-closed particle-hole 
lines (e.g. Fig.2a) and corresponding exchange terms Xi,X2, ■ ■ ■ (e.g. Figs. 3a, 5a), then: 



(i) it follows n{k) from the contraction SR (2.30) and from n{k) it follows t via (1.2). 

(ii) F{q) and P{r), which are needed for whe and gp{r), respectively, follow also from n{k), 
cf App.B. 

(iii) hg_p{r) and Ca,p(q') and the cumulant interaction energy Vc follow from Xa.,p- 

(iv) Finally, g^{r) = 1 - h^{q) and gp{r) = 1 - p{r) - hp{r). 

The resulting quantities are the GS-energies t, v = vup + Vc, e = t + v a.s functions of 
and the MD n{k) as well as the PDs ga,p{r) for the many-body quantum kinematics on the 
lowest level. 
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III. THE LOWEST-ORDER CUMULANT 2-MATRIX IN RPA 

In the following, these basic definitions and relations are applied to the spin-unpolarized 
HEG in the high-density limit — )■ 0, as an illustrative example. All linked-diagram (direct 
and corresponding exchange) terms are taken into account, which contribute in the energy 
to terms ~ r^pnrs + const + 0(rs)]. These diagrams are shown in Figs. 2a, 3a, 4a, 5a. 



3.1 The direct building block Xd in RPA 

Xdr denotes the direct cumulant 2-matrix Xd in its ring-diagram approximation of Fig. 2a: 
Xd = Xdr + ^Xd with an unkown (hopefully small) correction Axd beyond RPA. In its 
position representation, this Xdr is given by 

{/ S ^(?,^)<5i(g,r/)Q2(g,-r/)e^(^^^"'-^^^-'+9^-) + h.c.| 

(3.1) 

ri — r2 and Vur = ri — r\. The integral operator Qi{q, if) general- 



XdAri\r^,r2\r2) 



d^q 1 
47 2 



with the short-hands r 
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izes the particle-hole propagator Q{q, rj), cf (C.4). The diagonal elements Xdr(^i|^i7 '^21^2) 
^dr('"i2) define a function 



d^q sin qr 
47r qr 



du 



n 



R^{q,u)w{q,u) 



(3.2) 



which contributes with the RPA-beyond term A/ij to the CPD, = /i^r + A/id- When 



going from (3.1) to (3.2), the contour integration r] — )■ iqu, which replaces simultaneously 
the frequency integration by a velocity integration, leads to the real function R{q, u) = 



Q{q,iqu) for the particle-hole propagator, cf eg p2j, Eq.(B.l). The effective interaction 
w{q,u) = v{q,iqu) = v{q)/[l + v{q)R{q,u)] = q^/iq^ + qc^il^'^)] is the RPA screened 
Coulomb interaction with a Yukawa-like cut-off q'^R{q,u). The qualitative behavior of = 
hd is in [5j, Fig.3; it starts at /ia(0) = Ci and decays, passes a zero-intercept and approaches 



from below, its normalization is zero, see (B.l). The FT of /idr('"); namely 

00 

d^r 



Cdr{q) 



2tt 



q du R (g, u)w{q, u) 



(3.3) 



contributes to the CSF Cd = Cdr + ACd [29j with a beyond-RPA term ACd. How does 
Cdriq) behave for small and large momentum transfers, q qc and g — 00, respectively. 
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and in between ? 

For small q <^ gc, the dynamic RPA-propagator R{q,u) can be approximated by its static 
truncation Ro{u) = 1 — Marctan(l/M) |,29j such that Cdv{q ^ Qc) = L{y) q^jl with y = q/qc 
and 

oo 

y + -Ro{u) 

(For the properties of L{y) see [6].) Thus 

C,.iq « q.) = 4| - - f- V ■ ■ • ' (3.5) 

8 4wpi 4a;pi V?c/ 

in agreement with [35VI37] . The ellipsis represents terms originating from the difference 

AR{q,u) = R{q,u) — Ro{u). The function Cdriq) starts linearly in a small region (g ^ qc), 

which shrinks and finally vanishes for — )■ 0. The beyond- RPA factor Zp is added by hand 

with the following argument. As shown in the small-g behavior of the HF-function 



F{q) is given in (B.18). It transfers via the SFs to the CSFs and [together with the 
assumption ACd(g <^ qc) = —Zpq^/32 + 0(g^)] makes them behave correctly for small q, 
according to the plasmon SR (B.3) and ( B.13[ ), see also Table 1 and 2. 



At the other end, the large-g asymptotics does not require the ring diagram summation, 
so perturbation theory holds and the " descreening" replacement w{q,u) ~ v{q) is possible. 

oo 

With the Macke function /(g) = Svrg / du R'^{q,u) (explicitly given in p_j) it is 





"In between" g — )■ and g — )■ oo and according to (1.9), the integral Cdr = / d{q^) Cdr(g) 

I I 

contributes to the correlation parameter ci = Cdr + ■ ■ ■ , where the ellpipsis represents the 

missing terms beyond RPA. The behavior of Cd = Cdr + ^Cd with the RPA-beyond term 

ACd (causing the correlation parameters 01^2) is summarized in line 1 of Table 1 with 

unknown coefficients 54^5. 

3.2 The exchange building block Xx in RPA 

As stressed by Geldart et al. [33] , for a consistent small-r^ description the leading exchange 
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terms of n^{k), Sx{q) are needed. They follow from Xx with the ring-diagram part Xxr (Fig. 
3a): Xx = Xxr + ^Xx with an unkown (hopefully small) correction beyond RPA. That 
exchange contribution Xxrj which corresponds to the above Xdr, arises from the exchange 
r[ r'a, i.e. from Xxr{ri\r[,r2\r2) = Xdr{ri\r2,r2\r[), thus 

Xxrir,\r[,r2\r',) = ^ / ^ { / ^ Qiiq,v)Q2iq, -v) viq,v) e'(fc^^-'-^^^-'+9^-) + h.c.} 

(3.7) 

The diagonal elements h^ri^u) = Xxr(^i|^i; ^2|^2) define a function 

9 f d^Q f dfj ^ sin kv 

^"'^""^ " 2 y 17 2^^ ^i(^'^)^2(g,-r?) <q.v)^^\k=k,+k,+q ' (3-8) 

which contributes to the CPD hy^ = /ixr + A/ix. Its FT 

Cxr(g) = y ^ \h^r{r) = ^ y ^ Qi{k,v)Q2{k,-v) v{k^Tl)\k=k^+k2+q 

contributes to the CSF = Cxr + ACx with the RPA-beyond term ACx. A careful study 
shows, that Cxr can be written as 



^Pl f J3,,3, fkdu 6162 1 



^""■^^^ (47r/3)2 y ^ ^2 y ^ (^2^2+^2)^p^2 + ^2) " ^ + g2^(fc, «) ' ^^'^^ 

/ci,2 < 1, |fci_2 + g| > 1, fc = fci + fca + g, £1 = t(fc2 + q) - t{ki), £2 = i(fci + q) - t{k2) ■ 



This is the exchange counterpart to the direct CSF (3.3), unfortunately not known as an 
explicit function of q. If one assumes for the small-g behavior C^^iq ^ qc) = 0{q'^) and 
ACx(g ^ qc) = —c + 0(g''), then the correct small-g behavior of Cp = Cd — and Sp 



results, namely Cp(0) = c, see (B.2), and consequently 5'p(0) = 0, see (1.4). Line 2 of Table 
1 contains so far unknown coefficients S45. At the other end, the large-g asymptotics of 
Cxriq) can be extracted from perturbation theory, replacing "r" by its lowest order "1". 
This allows the w-integration, by which appears the well-known energy denominator: 

dx 6162 11 1 



vr (x2 + e2)(a;2 _^ ^2^ ei + £2 n + r2 q{ki + k2 + q)' 
Thus, it results similar as Cdv{q — ^ 00), namely 



CUq ^ 00) = + 0(c.p,), (4,/3^2 = ^ + ^ + ^ ) (3.10) 



with an unknown coefficient x. For I^iq) in general and possibly x = 8/5, see App.C. The 
behavior of = Cxr + ACx is summarized in line 2 of Table 1. An additional correlation 
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parameter arises from the RPA-beyond term ACx. 



3.3 How the cumulant structure factors Cap follow from C, 



For the spin-antiparallel structure factors it holds C-^ = Cd- In RPA, Cd{q) is given by (3.3) 
and its large-g asymptotics is given by (3.6). But beyond RPA the electron-electron coa- 
lescing cusp theorem h'^{0) = —ars[l — /ia(0)] (App.Bl, [31]) makes the first term decorated 
with an additional factor 1 — ci = 1 — hg,{0) < 1. One may expect the next term is modified 
in a similar way, thus 



00) = (l-Ci)^+(-+C2 



(3.11) 



with unknown correlation parameters C\ 2 vanishing for — )■ 0. Note that the spin- 
antiparallel CSF Ca(g — !■ oo) starts with 1/g'*. 



For the spin-parallel structure factors it holds Cp = Cd — Cx, Cpr = Cdr — Cxr, and ACp = 
ACd — ACx. The large-g asymptotics of Cx(g) is conjectured as (with another correlation 
parameter Cg) 

2 2 

Cx(g ^ oo) = (1 - ci) + (x + 4) ^ + ■ ■ ■ , (3.12) 
such that the difference Cd — Cx = Cp starts asymptotically according to 



Cp(g CX)) 



X + C2 - C2 



pi 



+ 



(3.13) 



with 1/g^ (not with ^lq^\ Just this conclusion results also from the electron-electron co- 
alescing curvature theorem with x = 8/5, cf ( |B.7 ). At the other end, the small-g behav- 
ior Cp(g ^ qc) follows from the above derived/assumed expressions for Cd{q ^ qc) and 
Cx(g < qc): 

4 / \ 5 



Cp — Cd — Cx — c -\- 



3q q 

y ~ 16 



q 



+ 



pi 



_q_ 



_q_ 
Wpi 



+ ■ ■ ■ 



c^pi. (3.14) 



[For the qualitative behavior of the corresponding CPD hp see [S], Fig. 4; it starts with 

hp{r < rj = (c3/2)r2 H , see ( |R2| ), ( [R5| , ( |R8| ).] A consequence for the total CSF 

C = Ca + Cp = 2Cd - Cx is 



C{q < gc 



^ 2 /^3g g3 

C + Zp 

^ V 4 16 



2uj 



+ 



pi 



C^pl 



+ 3 



Wpi 



+ ■ ■ ■ 



c^pi + ' 



(3.15) 
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This, used in S{q) = 1 - \F{q) - C{q) = c + 4(f - f^) - C{q) + cancels the 
first three terms. Thus emerges the quadratic term in agreement with the plasmon SR 
S{q <^ qc) = g^/(2cjpi) + ■ ■ ■ . Consequently an inflexion point appears near the origin. Its 
trajectory with — t- is sketched in App.B. Table 2 gives the small-g and large-g behavior 



of ^a^p . For the small-g behavior of S = Sg. + Sp the result (B.13) of the inflexion-point 



analysis with further correlation parameters 04^5 is used, see Table 2. It remains open to 
check the above assumption on Cx(g ^ qc)- 



3.4 How V follows from the cumulant structure factor C = Ca + Cp = 2Cd — Cx 



If (3.3) and (3.9) are used in (2.29), then the cumulant interaction energy Vc = Vd + fx 



^^dr + fxr H with ("r" = RPA, ellipsis = beyond RPA) 

1 f d?q 1 

Vat + v^r = j v{q) [ Cdr{q) - ^C^riq)] (3.16) 

consists of a "direct" part v^r and a corresponding "exchange" part v^r- The flrst one, 

00 

Vdr = -— I — q du v{q)R^{q,u)w{q,u) , (3.17) 



47r / 47r 







agrees (of course) with what follows from the total energy e after Macke (1950) [13] and 
Gell-Mann/Brueckner (1957) [I5], 



3 r r 

Cdr = / — g / du {v (q) R{q, u) -ln[l + v {q)R{q,u)]} 




(3.18) 



by means of the virial theorem v = Tsde/drg. The "exchange" part fxr can be simplifled as 



1 f d^q 



v{q)CAq) + 



-/x(g)f(g) . 



(3.19) 



A J An ■ ' (47r)2 

Namely, whereas in Cxr(q') the screened interaction v{q,r]) is necessary to remove the long- 
range divergencies of the bare interaction v{q), this is not the case for fx2, which has been 
calculated by Onsager/Mittag/Stephen (1966) ^8] with the result 



3 r d^q 

I J Wf 



^x2 = T / 7^ v\q)h{q) = 2 



6 4 7r2 



(3.20) 



Note that Vx2 = 2 6x2- 
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The agreement of (3.18) with the results of Macke and Gell-Mann/Brueckner [TU [15] and 



of (3.20) with the calculations of Onsager et al. |3H] confirms also the above expressions for 
hii,p{r) and Ca,p(g) in the RPA approximations "dr" and "xr". 



3.5 How follows from the contraction of Xx 



From XxT of (3.7) follows also a contribution to 

Xx(^) = Xxr(^) + ■ ■ ■ , Xxr(^) = 



r2 



Xxriri\r[,r2\r2) . (3.21) 



2 -3712 J 2-37r2 

As explained in App.D, the MD n{k) starts with a 2nd-order term, hence v{q,r]) can be 
replaced by —v{q)Q{q,ri)v{q,ri). Therefore 

9 f d^q 1 ( f dr] 



Xxr[ri\r^,r2\r2) 



4tt 2 



27ri 



7 [-v{q)Q{q,v)'v{q,v)] 



(3.22) 



■ Qi(g,^)Q2(g,-r?)e-'(^i+^^+^)^^ e^K^^+^^^^+^^^'^i +h.c.} . 

The r2-integration and afterwards the rnz-integration yield (including prefactors) 
6{ki + k2 + q) 6{k + k2)= i—j 6{ki + q-k)6{k2 + k) , 



2 ■ 3iT' 



so that 



1 r d^q 1 

2 / 17 2 



drj 
27ri 



v{q)Q{q, v)v{q, ri)X{q, r]) + c.c. 



with 



X{q,r])= Atx Qi{q,7])5{ki + q~k) ■ A-n Q2{q, -r])5{k2 + k) 



eik,q) 



(3.23) 



(3.24) 



(3.25) 



[T{k,q) -r]]"^ 

where 6(fc, g) = 1 for |fc| ^ 1, |fc+qr| ^ 1 and otherwise. The last step is derived in App.C4. 



This used in (3.24) gives for the MD-correlation part the result An{k) = ny.{k) + ■ ■ ■ with 
1 f d^q I 



n,{k ^ 1) = ±- 



^-(g)g(g,r/).(g,r,)^^^^^^y3^ 



+ c.c. 



2 J 47r 2 

With 0(fc, q) = ±0(fc, q) for |fc| ^ 1 and again with rj — iqu, it turns out 

oo oo +1 

nrik) = \ jd{q^) I g viq)R{q,u)w{q,u) Re J f p^^^^ 
-1 



(3.26) 



(3.27) 



being ^ for A; ^ 1. This is exactly the RPA-MD of Daniel/Vosko (1960) and Kulik (1961), 
as it has been shown in detail in p2] on the way therein from Eqs.(3.34) to (3.38). Also 
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therein at (3.10) it is shown, what rij-lk) contributes to the kinetic energy t in agreement with 



the virial theorem (1.1) and the Macke/Gell-Mann/Brueckner energy (3.18). The asymptotic 
behavior is n{k oo) = Loy{2k^) H for RPA and n{k ^ oo) = (1 - Ci)uj^^/ {2k^) H 



beyond it. This 2nd-order resuh is reduced by a factor of 1/2 through a still missing ladder 
diagram in its exchange version, as treated in the next section. 



IV. THE 2ND-ORDER 2-MATRIX IN RPA 

This ladder diagram has a particle-hole line running from 1' to 2 and a 2nd line from 2' to 
1, cf Fig.Sa. Its contraction r'2 = r2 and J (Pr2 yields the diagram of Fig.Sb, which differs 



from Fig.Sc only by a factor sign(l — k) according to (2.28). In detail the diagram Xha of 



Fig.Sa is given by (C.12). If for simplicity the RPA screening is neglected (corresponding to 



the replacement r— ;'2) and a series of laborious contour integrations are carefully performed 



together with the contraction steps analog (3.21), it finally yields 

d^qid^q2 ^pi ©(^,91,2 



> 



rv 



n^{k ^ 1) = TXix2(^) • 



(4.1) 



(47r/3)2 qfql {q, ■ q,y 

The mentioned factor sign(l — k) transforms Xix2(^) ^^e corresponding MD contribution 
rixik), which thus reduces both the particle branch {k > 1) and the hole branch {k < 1) of 



the RPA-MD nj.{k), Eq.(3.27). In particular the asymptotics for k ^ 00 changes from 1/2 
to 1/4 of ojpi/k^ as already mentioned at the end of Sec. III. Furthermore it is easy to show 
by means of the substitutions k — t- — (fc + q-^^) and ^2 ~Q2j (0 ^^ro normalization 
J (Pk nix2{k) = and (ii) the kinetic energy 



Ux2 



(47r/3)= 



(4.2) 



where also the identity (C.20) is used. The result (4.2) is just the exchange integral, which 



has been solved analytically by Onsager et al. [38]- Whereas (4.2) is a non-divergent 



2nd-order term, for n{k ^ 1) an indeed complicated ring-diagram summation is needed, to 
ameliorate diverging terms. 
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V. SUMMARY AND OUTLOOK 



This is a contribution to the theory of extended many-electron systems in terms of 
reduced-density matrices (RDMs), which replace $(1, • • • ,N), the many-electron wave 
functions of finite systems and form an infinite hierarchy, such that a lower-level RDM 
follows from the next-higher-level RDM by means of a certain contraction procedure, e.g. 
from the 2-body matrix (2-matrix) 72 follows the 1-matrix 71. This paper is driven by the 
behef that the cumulant 2-matrix 7c (dimensionless: x) has a fundamental meaning. 
This shows up for extended systems (i) in its size-extensive normalization and contraction 
and (ii) in its (closely related) /mA;eo?-diagram representation. Besides it is a decisive key 
quantity, because from it follows both the pair density (PD) and the 1-matrix. Here it 
is aimed to unveil its secrets for the high-density spin-unpolarized homogeneous electron 
gas (HEG). Although this model is only a marginal corner in the complex field of electron 
correlation, one should get generally deeper insight into the many-body correlation by 
exploiting the concept of cumulant decomposition of RDMs. The ground state (GS) energy 
of an electronic system is a simple functional of p{r, a) — electron density, of n{k, s) — 
momentum distribution (MD), and of 5(1,2) = pair density (PD) with 1 = (ri,(Ti). These 
1- and 2-body quantities have their common origin in the 2-matrix 72(1|1', 2|2'). Whereas 
the PD simply follows from 72 by taking the diagonal elements, the 1-body quantities need 
a procedure called contraction which means 2' = 2 and J d2. But this procedure leads in 
the case of extended systems to difficulties with the thermodynamic limit (TDL). Namely 
non-size-extensive expressions ~ appear, resulting from the numbers of pairs N{N — 1). 
This is avoided by the cumulant decomposition of the 2-matrix 72. This means that 72 is 
additively decomposed into a HF-like (or exact exchange) term 7hf, which reduces to a sum 
of products of the (correlated) 1-matrix, and a non-reducible remainder, called cumulant 
matrix 7c. Its advantage: it is size- extensively normalizablc and contractablc. Thus the 
problems with the TDL are eliminated. This is worked out in detail for the HEG with the 
homogeneous electron density p, with the MD n{k) and with S{q), the structure factor (SF) 
[or equivalently the 1-matrix /(|ri — r[\) and the PD g{ri2), respectively]. The 2-body 
quantities S{q) and g{r) are decomposed into HP terms and remainders C{q) and h{r), 
respectively called cumulant SF (CSF) and cumulant PD (CPD). In detail, this means 
g{r) = 1 — ^P{r) — h{r) and S{q) = 1 — \F{q) — C{q) with F{q)= Fourier transform of 
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p{r). The contraction procedure is developed in Eqs.(2.25)-(2.28): A quadratic equation 
for An{k) = n[k) — no(/c) yields the 2 branches for particles [k > 1) and holes {k < 1). 



Closely related to the described size-extensivity of the cumulant 2- matrix is its form 
(within many-body perturbation theory) in terms of linked diagrams. Under the impression 
that the triangle of "cumulants -H- size-extensivity -H- linked diagrams" is a fundamental 
relation (with an exponential-linked-diagram theorem as a more general mathematical 
background [101 H]), one may curiously ask for the cumulant geminals ^^(ri,r2) and 



their weights z/^, which diagonalize 7c |12], analogous to Eq.(2.16) in What are their 
properties? Is there perhaps a 2-body scheme, which allows to calculate approximately 
these ip"^ and z/^? How they differ from the Overhauser geminals [18], which parametrize 
the PD g{r)7 One access to this subject are the contracted Schrodinger equations [H]. But 
instead of this, in this paper the cumulant 2-matrix of the high-density HEG is derived. 
The results for 7c (dimensionless x) prove to be correct, because they yield the known 
RPA expressions for S{q), g{r), and n{k). This is the 1st step for further curiousity 
studies along the line of Overhauser [18] , namely to find the cumulant geminals of the high- 



density HEG, to be specific: the expressions (3.1), (3.7), and (C.12) have to be diagonalized. 



Furthermore, the coalescing cusp and curvature theorems, the plasmon sum rule (and its 
consequence, the inflexion-point trajectory), as well as the small- and large-g behavior of the 
CSFs for electron pairs with parallel and antiparallel spins are systematically summarized 
(analog to [36]), see Table 1 and 2. Within this, also the so far assumed asymptotic behavior 
of the structure factors for small and large q is corrected. 
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Appendix A: Generating functionals and linked diagrams 

The most compact quantity containing all the secrets of the HEG GS ) is the generating 
functional 

with X = {r,(j) and j dx = J d^f'- vi^) i]^{x) are Grassmann variables and P makes 
the normal ordering, i.e. orders the creation operators ip"^ including sign changes to the left. 
The Volterra coefficients of 7 [77^ 1 77] are the RDMs 71, 72, ... . The ansatz 7[^7^|77] = e^'^^'^^ 
defines another generating functional, xi'?^!^]? the Volterra coefficients of which are the 
so-called cumulant matrices xi, X2, ■■■ , such that a hierarchy results with higher level RDMs 
expressed by sums of products of lower level RDMs and non-reducible remainders. Thus 
it holds e.g. 7i(x|x') = 72(a;i|x'^, X2IX2) = Axi{xi\x[)xi{x2\x'2) + X2{xi\x[,X2\x2), 

with A being an antisymmetrizer [10]. Within time-dependent perturbation the- 
ory (S-matrix theory) the cumulant matrices Xi? X2, ■■■ are given by non-vacuum 
linked diagrams only and closely related to this they are size-extensively normalized i.e. 
J dxidx2 ■ ■ ■ xl^^il^^i) 3;2|a;2, ■ ■ ■ ) ~ A^, what allows for extended systems the thermodynamic 
limit with N ^ 00, Vt ^ 00, p = N/Q = const [41] . 



Appendix B: Structure Factors 

Here are summarized the normalizations of Ca,p(g), C{q) etc. as well as their behavior for 
small and large arguments. 

Bl: Perfect Screening Sum Rule. Cusp and Curvature Theorems 

S{q) starts with 5(0) = and approaches 1 for g — t- 00. g{r) starts with g{0) < 1/2 
and oscillatory approaches 1 for r — )■ 00. Note, that ga,^p{r) > rx g{r) > 0, since they 



are probabilities. The normalizations of 1 — S{q) and 1 — g{r) are contained in (1.5) with 
g{0) = (1 — ci)/2 (defining the correlation parameter ci) and with S{0) = (the perfect 
screening SR), respectively: 



00 



J d{q') [1 - S{q)] = 1 + ci , a' J d{r') [1 - g{r)] = 1 . 
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Similarly, the normalizations of 1 — ga.,p{r), 5'a(g), and 1 — Sp{q) are in (1.4) fixed with 

5'a(0) = Ci, 5'p(0) = 0, and 5'a,p(0) = 0. In terms of cumulant and "a,p" components, this 
reads as 

oo oo 

hM =ci^ j d{q^) CM = ci , CM = 0^a^ J d{r')hM = , (B.l) 



oo oo 

/ip(0) = O y" d{q') Cp{q) = , Cp(0) = c ^ a' J d{r^)h^{r) = c , (B.2) 



which follows from the FTs (1.9) and (1.10) for r = 0, respectively for q = 0. With 
|F(0) = 1 - c, the perfect screening SR ^(0) = is hidden in C(0) = ^(0) + Cp(0) = c. 
More general is the plasmon SR pTf 118] 

,2 



S{q < ge) 



2uj 



pi 



0{q' 



(B.3) 



see Table 2, last line. 



According to the electron-electron coalescing cusp theorem (7^(0) = avg 5'a(0) for the spin- 
antiparallel CPD [31], it holds h'M = -C(rs{l - ci) with 



ci = hM = -7- (tt^ -3 + 61n2) 
ovr 



2 3 



3« 
2^ 



r^lnr, + 0(r2) (B.4) 



Exactly this high-density behavior of h{fi) and g{fi) results also from the ladder 
theory as the best method to treat short-range correlation |l3l HH El]. By means of the 



Kimball trick pT] it follows the 1st term of Eq.(3.11), whereas the 2nd term follows from 



C&{q) = Cd{q) ~ Cdriq) and from (3.6). C2 is another correlation parameter vanishing for 
— )■ 0. Ci appears also in the large-Zc asymptotics n{k — )■ oo) = (1 — Ci)Upi/ (ik^) + ■ ■ ■ . 
Thus short-range correlation determine the large-wave number asymptotics of n{k) 
and of C{q) 



For the spin-parallel quantities, the Pauli principle makes 5'p(0) = 5'p(0) = or equivalently 
hp{0) = /ip(0) = 0. Furthermore the electron-electron coalescing curvature theorem says 
g'^'iO) = far, ^^^(O) or equivalently h'^'{0) = far, [2/"(0) + /i'^(O)]. With f"{0) = -ft and 

with m] 

C3 = K{0) = (47r2 - 5 + 201n2) ar, + ©(r^ Inr,) , (B.5) 
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it holds /ip'(O) = Q;rs(|c3 — 2t). This determines the large-g asymptotics of Cp{q) as shown 
in the following by means of the Kimball trick [21]. First a function Cp(g) ~ is defined 



through Cp{q) = a/{l + q^Y + Cp{q). Then from the FT (1.10) it follows 



/Zp(r) = a^(l + r)e-+ / d{q') c,{qf-^ 
16 J qr 



With the above curvature theorem it finally results 



h'i;{0) = a 



3tt 



(B.6) 



Cp(g^cx)) = (3 C3-4t)-f + 



(B.7) 



Comparing this with (3.13), the SR | — x + C2 — = Sc^ — At seems to hold. For — )■ it 



takes t — 3/10, hence x = 8/5 [for plausible arguments see text after (C.8) 



Because Cp{q — )■ oo) decays so strong, from (B.6) it follows in addition to (B.3) the relation 



C3 = h;{o) 



oo 

J dq q' CM) ■ 



(B.8) 



Since /ip(0) > 0, cf j32|,|17], the integral has to be negative. This is in agreement with (B.3) 



and (B.7): a factor q^ in front of Cp(g) makes the integral vanishing, cf (B.3). An additional 



factor q enhances the negative asymptotic branch of Cp(g — )■ oo), cf (B.7). Besides, from 
gl{Q) > it follows /i'j^(O) < 4t/3. For t(r,) cf eg [22]. 



B2: Plasmon SR and inflexion-point trajectory 



In (3.5) it has been shown how perturbation theory and RPA makes the plasmon-term ~ q 



arise from the linear term ~ g of the non-interacting system. Here it is asked on the contrary, 
if the interaction is switched off (r^ — )■ 0), how the quadratic behavior ~ q^ is transformed 
to the linear behavior ~ q. The answer: if S{q) starts for small q according to the plasmon 
SR [18] quadratically with g^, then there must be - already from a topological point of view 
- an inflexion point go) 'S'(go); which moves with — )■ towards the origin. So also the 
linear behavior at the inflexion point is transported to the origin, realizing thus the correct 
behavior of the unperturbed ideal Fermi gas, as it should. Now this qualitative discussion 
is quantized by the ansatz 

-,2 



S{q < ge) 



2uj 



pi 



-[1 - az^ + hz-^ + 



q_ 



(B.9) 
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and adjusting the r^-dependent coefficients a, b appropriately. Namely the power expansion 
of S{q) around g ~ go 

S{qo + x) = S{qo) + ^S\qo)x + ^S"{qo)x^ + ^S"'{qo)x' + ■■■ (B.IO) 

is compared with its "r^ — )■ 0" limit So{x). This gives 3 equations, which allows determining 
the 3 quantities go, cl, b: 



^S\qo) ^ , ^S"{qo) = , ^S"'{qo) ^ 



with the result (for — > 0) 



Wpi 2 



1 , 22 
a — )■ — , — 7- — — 
32 ' 33-5 



(B.ll) 



(B.12) 



So (B.9) takes for Ts — ?■ the form (with unknown correlation parameters 04,5) [35] 

2 



^(g < gc 



2w 



pi 



1 



1 



C4 



+ 



(B.13) 



Thus the inflexion point moves towards the origin according to go — )■ 1.5 Up\ and S'(go) — >■ 



0.96 cjpi with a slope of 0.64 being smaller than S'g(O) = 0.75. Comparison of (B.13) with 



(3.15) shows S4" = ^(^- + 04) and s'^' = "^(pTs+Cs). - As quoted in [33], other consequences 



(B.14) 



of (B.9) are g{r > Tc) = 1 - 2V(c^pi r^) ^ and 

00 00 

J dr /[I - g{r)] = — , J dr r^[l - g{r)] 


The equation with the weight is widely used in [37]. Remind 1 — g{r) = ^/^(r) + h{r), 
what relates J d^k {dn{k)/dky to J d^r r'^h{r) for 1/ = 2 and 1/ = 4. 
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B3: The HF-function F(g) 

For the spin-parallel quantities gp{r) ^H- 5'p(g) and hp{r) ^ C'p(g) the functions P{r) and 



its FT F{q) are needed, cf (1.7), (1.8) and Table 2: 



F{q) = a' [d{r')'-^f{r) , 1^(0) = 1 - c , f d{q^) If {q) = I 



qr 


F(g) transports the MD n(/c), as seen from 

1 . 3 



(B.15) 



2^(^) = ^ / d'kn{k)n{\k-q\) 



(B.16) 
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Thus it is related to the probabihty of finding a pair of electrons with given relative mo- 
mentum q |3n]- For = with no(fc) = 0(1 — k) the integral takes a simple geometric 
meaning: it is just the volume of two spherical calottes with the height h = 1 — g/2, thus 
/ d3A;... = 2-f/i2(3-/i) or 

oo 

0(2 -g), I d{q') ^-Fo{q) = 1 . (B.17) 



[By the way, (i) -Fo(g) appears in the Overhauser theory as the weight of its PD geminals 
[15] and (ii) the singularities of F{q) at g — and g — > 2 are studied in [S3]-] The difference 
AF(g) = -F(g) — FqIq) results from An{k) = n{k) — nQ{k). It holds 

^AF(g«g.) = -c-4(f -f^)+ ^AF(g oo) = + • • • , (B.18) 

oo 

and / d{q^)AF{q) = 0. 



B4: The jump discontinuity of S"{q) and C"{q) at g = 2 

For transition momenta \q\ = 2, when passing this value from below, the topology changes 
from 2 overlapping to 2 non-overlapping Fermi spheres. The consequence: A/"(2) = 27r^, cf 
(C.2) in [9J r\ AC"(2) = (3ci;pi/4)^ + ■ ■ ■ [6J. The ellipsis represents exchange and beyond- 
RPA terms. For S"'(2) the relation lAF"(2) = -3z|/4, cf (5.2) in [53], has to be taken into 

account. Thus in lowest order AS'"(2) = 32:^/4 — (3(X'pi/4)^ H , what causes (in the direct 

space) the Friedel oscillations of g{r). It remains open to study the exchange term A/"(2), 
which presumably reduces the direct term A/" (2). 



3g 1 /g^3 
22 2 V2. 



Appendix C: Particle-hole propagator and Macke function I{q) 

Cl: The building elements of the RPA Feynman diagrams (cf Figs. 2a, 3a, 4a, 5a) are the 
Coulomb repulsion t>(g) = g^/g^ with the coupling constant g^ = iars/n and the one-body 
Green's function of free electrons with t{k) = k'^/2, 

^ , 0(A;-1) 0(1-A;) 

Go{k,uj) = \ ' + \ ' , 5>0. Cl 

u — t[k) + id u — t[k) — lb 

C2: From Go{k,uj) follows the particle-hole propagator Q{q,ri) in RPA according to 

/H^k f dill 
— J —Go{k,uj)Go{\k + qluj + v) (C.2) 
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with the result 



47r 



e+(fc,g) 



7] — T{k, q) — i6 T] — T{k, q) + i5 
T{k, q) = t{k + q)- t{k) = g (^A;C + I) , C = cos <(/e, q) . 



(C.3) 



The denominators contain the excitation energy T{k,q), to create a hole with k inside the 
Fermi sphere and a particle with k+q outside the Fermi sphere. A complicated step function, 
defined by 6^ = 1 for k ^ l,|fc + qr| ^ 1 and otherwise, is a consequence of the Pauli 



principle. R{q,u) = Q{q,iqu) defines a real function. A generalization of (C.3) is 



Qi{Q,v)f{ki) 



An 



Q1 



e; 



r]-T,-i6 r]-Ti + i6 



f{ki) , Ti = T{ki,q) 



defining the integral operator Qi{q,ri), i = 1,2. 



(C.4) 



C3: When calculating the CSF Cid, then the Macke function I{q) appears via 



dr] 



d^kid^k2Qi 



(C.5) 



27ri ^ " (47r)2 ' J q . + k2 + q) 

When dealing with the corresponding exchange term Cix, then by means of contour inte- 
gration it results 



dr] - 2 - 

where the integral operator /(g) and the exchange term Ix{q) are defined by 

d^kid^k2Q)i2 



(C.6) 



I{q)v{ki + k2 + q) 



q-{ki + k2 + q) 

I{q) and I^iq) have the large-g asymptotics 

/(g — )■ oo) / 1 2 1 \ /x(g -> oo 



v{\ki + k2 + q\) = h{q)v{q) 



{An/sy 



1 2 1 



i X 
g^ g* 



(C.7) 



(C.8) 



Whether the guessed value of x = 8/5 [text after (B.7)] is correct, has to be studied. The 
change of 2/5 to 8/5 is a heavy enhancement of the 1/g+tail, which is plausible in view of 
the following qualitative arguments. In [9] and [12], the (explicitly not known) function 
/x(q') is compared with the explicitly known Macke function I{q). Whereas /(g) starts 
linearly with g (what causes the Heisenberg divergence), its exchange pendant /x('?) must 

oo 

start at least with g^. Otherwise the integral / dq Ix{q)/q'^ (which is proportional to the 
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congenially calculated Onsager-et-al energy Vx2) would not exist. An estimate shows even 
/x(g — > 0) ~ q^. This flattening in the small-g region together with the "area conservation" 

oo 

J dq [/(g) — /x(g)] = enforces the mentioned tail enhancement. 



C4: Proof of (3.25): In the following the notation 



e+(fc, q) = &{k - 1)0(1 -\k + q\), e-(fc, q) = 0(1 - A;)e(|A; + q\ - I) 



is used. (C.4) inserted in (3.25) gives 



X{q,v) 



eih - 1)0(1 - |fci + q\) 0(1 - A;i)0(|fci + q\-l 



T] - r(fci,qr) - i6 



T] - T{ki,q) +iS 



ki^k q 



(C.9) 



0(A;2- 1)0(1 -|fc2 + g|) _ 0(1 - A;2)0(|fc2 + g| - 1 

-V-r{k2,q) -i6 -7] - T{k2,q) + i6 jk2^_k 

0(|fc 1)0(1 -fc) _ 0(1- \k-q\)e{k-l ' 
7] -T{\k - q\,q) -i6 r] - T{\k - q\,q) + i6 _ 
Q{k - 1)Q{1 - \k - q\) _ Q{1 - k)Q{\k - q\ - 1 
—1] — T{k, —q) — \5 —1] — T{k, —q) + i5 

Only the products "over cross" contribute. With the substitution q — )■ —q it results 

+ I e(l-fc)0(|fc + g|-l) _ 0(fc,g) 



(C.IO) 



(C.ll) 



[r, + T{k,q)Y ' [r^ + T{k,q)Y [r/ + r(fc, g)]^ ' 

To have finally exactly the same expression as in [12] the substitution fc + g — > —k' changes 



the sign in the dominator, proving thus (3.25). 



C5: Proof of (4.1): The diagram of Fig.5a means in detail 



If d?qid^q2 f d^kid^k2 

i|e'l(fe+'Z,+g,)i-„,+fe,r-„,+(q,+q,)r„l x-,_2(fci,2, 91,2) + h.c.} 

where the following abbreviations are used: 

duiduj2driidri2 



xi{ki) 

X2{k2) 



Xlx2(fcl,2,gi,2) , ^27ri)4 

(fci^l) i\k, + q,\^l) 



(C.12) 



Xi(fci) ■ ^2(^2) with 

(1^1+^1+921 ^ 1) 



(C.13) 



LOi - t{ki) ±16 LOi+rji- t{ki + g J ±15 UJi + rji - r]2 - t{ki + qi + ^2) + i*^ 

{h ^ 1) (|fc2 + g2l ^ 1) (|fc2 + gl+g2l ^ 1) 

002 - t{k2) ±16 002 -r]2- t{k2 + ^2) ± i*^ 002 + 111 - r]2 - t{k2 + qi + ^2) + i'^ 
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For the contribution of the ladder diagram x\^2 to the MD n{k)^ the contraction steps analog 



(3.21) can be performed with the result 



^^^11' -ikr„, f ^if(fc2+q,+qjri2+fcir,,,+(q,+q.,)ri 
2-37r2 



2 ■ 37r2 



Air 

Y 



5(fcl-fc)5(fc2-fcl) 



Xix2(fc7 9i2) contains contour integrations in 4 complex frequency planes. Unfortunately, 
each of these 4 integrals runs over 3 frequency denominators in a difficult way. But by means 
of a trick ( integration variables are substituted) the ?7i 2 integrations can be decoupled: 
U2 0J2 ~ Tji + ri2. This does not influence Xi{k) (it is only written in reverse order), but it 
changes X2{k): 

(|fcl+ 91 + 921^1) + (fcl^l) 



Xi(fci) 
X2(k 



Ui + r]i - r]2 - t{ki + qi + 92) ±15 uji+rji- t{ki + q^) ±15 Wi - i(:(fci) ± \5 

(C.15) 



^ 1 



|fc + q2| ^ 1 



1^ + ^1 + 921 ^ 1 



1^2 - ^71 + ^72 - ^(fe) ± W U2- rii~ t{k + q2) ±15 U2 - t{k + qi + qs) ± i^ 

(C.16) 



Now it is easy to perform the integration over rj = t]2 — rji, namely 



(|fc + qi + q2| ^ 1) 



(A:^l) 



{T){k ^ l,|fc + qi + q2| ^ 1) 



drj 

27ri Ui — r] — t{k + qi + q2) ± i^i UJ2 + r] — t{k) ± 182 (jJi+ UI2 — t{k) — t{k + qi + q2) ± 1(5 

(C.17) 

Therein those cases do not contribute, where the 2 poles in the complex ?7-plane are both 
on the same side of the real axis. Only the combinations (±i5i) with (±1^2) contribute. 
Similarly the ?7i-integration is handled: 



{\k + q,\^l) 



(1^ + 921 ^1) 



{T){\k + q.2\^r 



dr]i 

2tx\ uji + Tji — t{k + q^) ± i5i UJ2 — rji — t{k + q2) ± 1^2 001+002 — t{k + q^) — t{k + q2) ± i5 

(C.18) 



The next step would be the Ci;2-integration. The integrand is a product of the rhs's of (C.17) 



and (C.18) with the 3rd factor of (C.16): 



duo2 _ {T){k^ l,|fc + qi + q2| ^ 1) 
27ri U01+ b02— t{k) — t{k + qi + q2) ± i^ 

(T)(A:^l,|fc + qi^2l^l) 
Wi + (X'2 — t{k) — t{k + q^) — t{k + q2) ± i5 

. (|fc + 91 + ^21^1) 

U2 - + ^1 + ^2) =t 



(C.19) 



0^(fc,9l,2) 



[uoi - t{k) + q^ • q2 ^ 15] q^ ■ qg 
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The integration region is defined by the functions 2) = 1 ^ 1) 1^ + 9i + 

^2! ^1,1^ + ^12! ^ 1) otherwise. The identity 



t{k) + t{k + q^ + q^) - t{k + q{) - t{k + q^) = qi ■ q2 



(C.20) 



is used. The remaining cji-integration finally yields 
dui Q^{k,q,2) 



Xlx2(fc,gi,2) 



Q(fc,gi,: 

27ri [ui - t{k) ± 16] [ui - t{k) + q-^-q^^ iS] q^ ■ q^ {q^ ■ ^2) 



(C.21) 



where 9 = 0+ + 0". (C.21) inserted in (C.14) gives (4.1), qed 



Appendix D: Correlation parameters 

Here is a list of characteristic correlation parameters as functions of r^, vanishing for — )■ 0: 

= d{k^) [n{k) — riQ^k)] is the normalization of the MD-correlation tails, 
c = dijc') n{k)[l — n{k)\ = N — d{k^) [n{k) — riQikyf' is the Lowdin parameter. 
1 — 2;f is the reduced discontinuity jump of n{k) at the Fermi surface |fc| = 1. 
s = — diJc") [n{k) \n.n{k) + (1 — n{k)) ln(l — n{k))\ is the particle-hole symmetric corre- 
lation entropy, see Eq.(22) in [22] . 

Ci = /ia(0) with 1 — /?-a(0) = S'alO) being the on-top value of the Coulomb hole, 
Ci simultaneously determines the large-g asymptotics of Cd,x(Q') and n(/c), 
C2,C2 determine the higher-order large-g asymptotics of Cd,x(q'), 

C3 = /ip(0) and t determine the on-top curvature of the Fermi hole 5'p(0) = |t — /ip(0), 
C4,5 determine the small-g behavior of S{q) beyond the plasmon-SR term g^/(2u;pi), 
Cg influences the jump discontinuity A5"'(2). 

All these parameters may serve as measures of the correlation strength. They vanish for 
— 7- with terms ~ r^, Inr^ or rl^"^ . From a diagramatical point of view all the quantities 
characterizing n{k) or derived from it, should start with r^. The reason: the diagrams for 
n{k) start with n2{k), the Ist-order term ni{k) (Fig. 4c) is zero - in agreement with the 



vanishing Ist-order kinetic energy ti = 0, what follows from the virial theorem (1.1). But 
there are starting dependencies as A = 0(rs^^) cf (3.46) in [12j. This peculiar behavior is 
an example for how higher-order partial summations may create lower-order terms. It is a 
peculiarity of the Coulomb repulsion. A ~ Vs^'^ results from that fc-region — 1| ^ Qc, where 
the RPA reconstruction takes place. A characteristic non-analyticity is also the behavior of 
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n{k) at the Fermi surface |A;| = 1, namely ~ (fc — 1) In |fc — 1 



[1] M.P. Tosi in: N.H. March (Ed.), Electron Correlation in the Solid State, Imperial College 
Press, 1999, p. 1. - G. Giuliani and G. Vignale, Quantum Theory of the Electron Liquid, 
Cambridge University Press, Cambridge, 2005. 

[2] There is a non- vanishing probability of finding 0, 2, 3, ... electrons in such a Wigner sphere. 
Such particle number fluctuations in fragments have been studied by P. Ziesche, J. Tao, M. 
Seidl, and J. P. Perdew, Int. J. Quantum Chem. 77, 819 (2000) with the conclusion 'correlation 
suppresses fluctuations', cf also P. Fulde, Electron Correlations in Molecules and Solids, 3rd 
ed., Springer, Berlin, 1995, p. 157. 

[3] The Coulomb repulsion takes the form arg/r, if lengths and energies are measured in units 
of k^^ and kp, respectively. This shows, that is not only a density parameter [2], but 
also measures the interaction strength or coupling constant. The FT of ar^jr is ^■nargjc^ = 
7r^(g^/(7^) with q measured in units of kp and = 4ars/7r. - Notice that e = k^e = ejiarsf' 
gives the energy in a.u., e.g. the energy in zeroth order and the lowest-order exchange energy 
are eo = 3/(10 a^r^) and Cx = — 3/(47rars), respectively. 

[4] J. Cioslowski, P. Ziesche, and K. Pernal, Phys. Rev. B 63, 205105 (2001); J. Chem. Phys. 
115, 8725 (2001). 

[5] P. Ziesche, Int. J. Quantum Chem. 90, 342 (2002). The last sentence of the Abstract has to 
be deleted. 

[6] P. Ziesche and J. Cioslowski, Physica A 356, 598 (2005). - Here it is shown, which peculiarities 
of n(A;) and Sijq) caused by RPA lead to the logarithmic non-analyticity of the high-density 
correlation energy Ccorr s l^^s" ~ Eq. (27) th.6 correct prefcictor is '^OLVgji^'K^. 

[7] R. Diez Muiho, I. Nagy and P.M. Echenique, Phys. Rev. B 72, 075117 (2005). 

[8] J. Cioslowski and P. Ziesche, Phys. Rev. B 75, 085103 (2007). Here the self-energy S(l, 1/2) 
is investigated. 

[9] P. Ziesche, phys. stat. sol. (b) 244, 2022 (2007); the relation /i - = S(l,/x) has been 
studied for — )• 0. - In Appendix C the Macke function /(g) is explicitly given. The notation 
/(g) is related to the notation used in [29] as A5'(g/2)D = — 3/(87r^)(gc/(7)^/(g) and in [6] as 
K{q) = 3/{87rV)I{q)- 



33 



[10] P. Ziesche, Ann. Physik (Leipzig) 16(1), 45 (2007). Here it has been shown S2x(l, 1/2) = e2x- 
[11] M.L. Glasser and G. Lamb, J. Phys. A 40, 1215 (2007). 

[12] P. Ziesche, Ann. Phys. (Berhn) 522(10), 739 (2010). - In the 2nd hne of Eq. (C.19) there is a 
prefactor missing. Read the last line on p. 762 as "respectively. One may guess, that /x(9) 
starts with ~ g^." . 

[13] W. Heisenberg, Z. Naturf. 2a, 185 (1947). 

[14] W. Macke, Z. Naturf. 5a, 192 (1950). - Another access with the same result is due to D. Bohm 

and D. Pines, Phys. Rev. 92, 609 (1953) and D. Pines, Phys. Rev. 92, 626 (1953). 
[15] M. Gell-Mann and K. Brueckner, Phys. Rev. 106, 364 (1957). 

[16] N.H. March, Phys. Rev. 110, 604 (1958). - Note that the "Hellmann-Feynman theorem" had 
been formulated first by P. Giittinger, Z. Phys. 73, 169 (1932); later it was to be found in H. 
Hellmann, Einfiihrung in die Quantenchemie (Deuticke, Leipzig, 1937), pp. 61 and 285 as well 
as in R.P. Feynman, Ph.D. Thesis, MIT, 1939; Phys. Rev. 56, 340 (1939). 

[17] D. Pines and P. Nozieres, The Theory of Electron Liquids, Vol. I (Benjamin, New York, 1966). 

[18] From the /-SR or Thomas-Reich-Kuhn SR of the dynamic SF S{q, oj) follows the plasmon SR 
S{q < gc) = gV(2wpi) H of the static SF S{q). 

[19] E. Daniel and S.H. Vosko, Phys. Rev. 120, 2041 (1960). 

[20] I.O. Kulik, Zh. Eksp. Teor. Fiz. 40, 1343 (1961) [Sov. Phys. JETP 13, 946 (1961)]. 
[21] P. Ziesche, phys.stat.sol. (b) 232, 231 (2002). 

[22] P. Gori-Giorgi and P. Ziesche, Phys. Rev. B 66, 235116 (2002). - An anaZytzca/ expression for 
n{k) in the range Vg = 1, • • • , 10 is carefully fitted to the well-known kinetic energy t{rs). - 
For the momentum distributions of the spin-polarized electron gas see [23]. 

[23] P. Ziesche and F. Tasnadi, Ann. Physik (Leipzig) 13, 124 (2004), Erratum 13, 623 (2004), 
Sec. 4, Fig. 2. - Besides in Sec. 5 it is also presented the spin-structure of the spin-polarized 
HEG. 

[24] M. Holzmann, B. Bernu, C. Pierleoni, J. McMinis, D.M. Ceperley, V. Olevano, and L.D. Site, 

Phys. Rev. Lett. 107, 110402 (2011) with Monte Carlo calculations for n{k). 
[25] P. Ziesche, Thesis, TH Dresden, 1963; Ann. Phys. (Leipzig) 21, 80 (1968). 
[26] LA. Nechaev and E.V. Chulkov, Phys. Rev. B 73, 165112 (2006). 
[27] A.J. Glick and R.A. Ferrell, Ann. Phys. 11, 359 (1960). 
[28] D.J. Geldart, Can. J. Phys. 45, 3139 (1967). 



34 



[29] J.C. Kimball, Phys. Rev. B 14, 2371 (1976). 

[30] D.M. Ceperley and B.J. Alder, Phys. Rev. Lett. 45, 566 (1980), with Monte Carlo calculations 
for g{r) . Other papers are based eg on the dielectric function £{k, uj) and local field corrections. 
[31] J.C. Kimbah, Phys. Rev. A 7, 1648 (1973). 

[32] A.K. Rajagopal, J.C. Kimbah, M. Banerjee, Phys. Rev. B 18, 2339 (1978). 

[33] D.J.W. Geldart, A. Houghton, S.H. Vosko, Can. J. Phys. 42, 1938 (1964). 

[34] P. Ziesche and F. Tasnadi, Ann. Phys. (Leipzig) 13, 232 (2004); Erratum 13, 624 (2004). 

[35] N. Iwamoto, Phys. Rev. A 33, 1940 (1986), Eq. (4.2). 

[36] P. Gori-Giorgi, F. Sacchetti, G.B. Bachelet, Physica A 280, 199 (2000), Table on p. 201; Phys. 
Rev. B 61, 7353 (2000); Phys. Rev. B 66, 159901(E)(2002). 

[37] L. Zecca, P. Gori-Giorgi, S. Moroni, and G.B. Bachelet, Phys. Rev. B 70, 205127 (2004). 

[38] L. Onsager, L. Mittag, and M. J. Stephen, Ann. Physik (Leipzig) 18, 71 (1966). Here the 
exchange integral of the 3-dimensional HEG has been evaluated analytically. This rather her- 
culean work was generalized and extended to the d-dimensional electron gas by M.L. Glasser 

m- 

[39] M.L. Glasser, J. Comp. Appl. Math. 10, 293 (1984). 
[40] P. Ziesche, Commun. Math. Phys. 5, 301 (1967). 

[41] P. Ziesche, in: J. Cioslowski (Ed.), Many-Electron Densities and Reduced Density Matrices, 

Kluwer Academic/Plenum Publishers, New York, 2000, p. 33. 
[42] P. Ziesche, Lecture Series on Computer and Computational Sciences 7b, 1051, (2006). 
[43] Z. Qian, Phys. Rev. B 73, 035106 (2006). 

[44] J. Cioslowski and P. Ziesche, Phys. Rev. B 71, 125105 (2005); Phys. Rev. B 72, 

239901(E)(2005). 
[45] J.C. Kimbah, J. Phys. A 8, 1513 (1975). 
[46] H. Yasuhara and Y. Kawazoe, Physica 85A, 416 (1976). 
[47] V.A. Rassolov, J.A. Pople, and M.A. Ratner, Phys. Rev. B 62, 2232 (2000). 
[48] A.W. Overhauser, Can. J. Phys. 73, 683 (1995). Succeeding papers are [igtiSH] . 
[49] P. Gori-Giorgi and J. P. Perdew, Phys. Rev. B 64, 155102 (2001): Extended Overhauser model. 

P. Gori-Giorgi and A. Savin, Phys. Rev. A 73, 032506 (2006): Electron-electron coalescence 

and short- and long-range correlation. 
[50] M. Corona, P. Gori-Giorgi, and J.P. Perdew, Phys. Rev. B 69, 045108 (2004). 



35 



[51] R. Asgari, M. Cardenas, M. Polini, B. Davoudi, and M.P. Tosi, |cond-mat/0408293 
[52] P. Ziesche, Phys. Rev. B 67, 233102 (2003). 

[53] P. Ziesche, phys. stat. sol. (b) 241, 3544 (2004): I; (b) 242, 2051 (2005): II. - It is shown, 
that the Kimball-Overhauser geminals do not satisfy the plasmon SR. 

[54] H. Yasuhara, Sohd State Commun. 11, 1481 (1972); Physica 78, 420 (1974); J.M. Herbert, 
Phys. Rev. A 66, 052502 (2002); R. Asgari, M. Pohni, B. Davoudi, and M.P. Tosi, Sohd State 
Commun. 125, 139 (2003). 

[55] L. Kong and E.F. Valeev, J. Chem. Phys. 134, 214109 (2011). 

[56] D.R. Alcoba, R.C. Bochicchio, L. Lain, and A. Torre, J. Chem. Phys. 133, 144104 (2010). 

[57] K.R. Shamasundar, J. Chem. Phys. 131, 174109 (2009). 

[58] S.D. Hamieh and H. Zaraket, Eur. Phys. J. D 55, 229 (2009). 

[59] M. Nooijen, M. Wladyslawski, and A. Hazra, J. Chem. Phys. 118, 4832 (2003). 

[60] D.A. Mazziotti (Ed.), Reduced Density Matrices: With Applications to Many-Electron Atoms 

and Molecules, Wiley-Blackweh, Hoboken, 2007; T. Juhasz and D.A. Mazziotti, J. Chem. 

Phys. 125, 174105 (2006). 
[61] J.M. Herbert, Int. J. Quantum Chem. 107, 703 (2007). 

[62] J.M. Herbert, "Reconstructive Approaches to One- and Two-Electron Density Matrix The- 
ory", Dissertation, Uni of Wisconsin-Madison, 2003, p. 88. 
[63] J.E. Harriman, Phys. Rev. A 65, 052507 (2002). 

[64] W. Kutzelnigg and D. Mukherjee, J. Chem Phys. 110, 2800 (1999); Chem. Phys. Lett. 317, 
567 (2000). 

[65] W.H. Young and N.H. March, Proc. Roy. Soc. A 256, 62 (1960). 

[66] M. Nakata, B.J. Braams, K. Fujisawa, M. Fukuda, J.K. Percus, M. Yamashita, and Zh. Zhao, 

J. Chem. Phys. 128, 164113 (2008). 
[67] J. Sun, J.P. Perdew, and M. Seidl, Phys. Rev. B 81, 085123 (2010). The correlation energy is 

interpolated between high- and low-density limits. 
[68] P.-F. Loos and P.W. GiU, Phys. Rev. B 84, 033103 (2011). Here, the spin-polarized HEG at 

high-density is treated. 

[69] The GS-energy e is also a functional of t(fc) and v{q) with 5e/5t{k) ~ n{k) fT9] and 6e/6v{q) ~ 

S{q) - 1 [291110]. 
[70] H. Yasuhara, Physica 78, 420 (1974). 



36 



Figures 
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FIG. 1: Feynman diagrams for the Coulomb repulsion in RPA with the partial summation (or 
screening replacement) v{q) — t- v{q, rj), q = momentum transfer, r] = energy/frequency, dashed line 



bare interaction v{q) = q^/q^, closed loop = particle-hole propagator Q{q,ri) of (C.3), wavy line 
effectively screened interaction v{q,r]) = v{q)/[l + v{q)Q{q,r])]. 
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FIG. 2: 2a = lowest-order renormalized cumulant 2-matrix XdD 2b = non-divergent cumulant 
PD hdr or cumulant SF Cdr (Kimball) |31] . 2c = non-divergent interaction energy f^r (Macke, 
Gell-Mann/Brueckner) [HI [15]. 
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FIG. 3: The RPA exchange terms corresponding to Fig. 2 with 3a = Xxr (following from 2a = Xdr 
through the exchange replacement 1' -H- 2'), 3b = /ixr or C^t (non-divergent), and 3c = t^xr- As 
shown by Onsager et al. |38j . already Vx2 does not diverge (needs no RPA renormalization) . 
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FIG. 4: How the 2-body diagram 4a = 3a = Xxr transforms to the 1-body diagram 4b through 



the contraction [with 2'=2 and / d2, see (2.25)], marked by a small circle and finally to the non- 



divergent RPA-diagram nj-{k) of Fig. 4c (Daniel/Vosko and Kulik) |19| l20j. The Ist-order term is 
zero (analog to the Ist-order kinetic energy ti = 0). 
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FIG. 5: How the 2-body ladder diagram 5a in its exchange version xixr transforms to the 1-body 



diagram 5b through the contraction (2.25) and finally to the diagram nx{k) of Fig. 5c. In Sec. IV 



and App.C5 the RPA screening is neglected for simplicity corresponding to the replacement r— )-2 
or v{q, rf) — v{q). 
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ha(r) 



■> ga(0 



h(r) ► V, 



ho(r) 



n(k) 



I / 



f(r) ► Vhf 







FIG. 6: Flow charts showing the way from the starting point Xd = direct cumulant 2-matrix 
through its exchange pendant Xx and the cumulant PD h(r) to the final results n(A;) and ^(r), 
from which follow (i) /(r) = 1-matrix, t = kinetic energy, c = Lowdin parameter and (ii) v 
= interaction energy, respectively. "a"=spin-antiparallcl, "p"=spin-parallel. Prom FT arises a 
similar scheme for the structure factor S{q), the cumulant structure factor C{q), and the HF term 
F{q) g{r), h{r), and f^{r), respectively]. Note that g{r) = ^[ga, + gp{r)] and e = t + v. 
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